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Temperature-dependent magnetization measurements on a series of synthetic ferritin proteins con- 
taining from 100 to 3000 Fe(III) ions are used to determine the uncompensated moment of these 
antiferromagnetic particles. The results are compared with recent theories of macroscopic quan- 
tum coherence which explicitly include the effect of this excess moment. The scaling of the excess 
moment with protein size is consistent with a simple model of finite size effects and sublattice 
noncompensation . 



PACS numbers: 73.40. Gk, 75.60.Jp, 75.10.Jm 

In nanometer-scale magnetic particles it is possible to 
observe a number of phenomena which do not exist in 
bulk samples. The evolution of magnetic order from in- 
dividual atoms to large clusters, Id thermal rjelaxation of 
the magnetization from a metastable statc,EI the differ- 
ent roles of surface and interior atoms, a and the|-jquan- 
tum mechanical dynamics of the order parameteio have 
all received considerable experimental and theoretical at- 
tention in the last several years. Here we present mea- 
surements on a series of well-characterized samples of 
biomimetic antiferromagnetic particles with sizes from 
100 to 3000 Fe(HI) ions per particle. By separating the 
bulk from the surface contributions to the magnetization, 
we explore the connection between two of these phenom- 
ena: the role of the excess moment in the macroscopic 
quantum coherence of antiferromagnetic particles. 

Past theoretical work suggested that the Neel vector of 
a small antiferromagnetic particle could exhibit macro- 
scopic quantum coherence (MQC), in which it tunnels 
resonantly between degenemte easy directions, at a rate 
accessible to experiment .ErQ Measurements of the mag- 
netic noise spectrum and ac susceptibility of the antifer- 
romagnetic cores of the protein ferritin revealed a reso- 
nance whose frequency scaled with particle size, applied 
magnetic field, temperature, and interparticle separatipu 
in qualitative agreement with theoretical predictions. lTeI 
In the interpretation of these results, it was assumed 
that any excess moment of the ferritin cores would fol- 
low the dynamics of the Neel vector without affecting it. 
Measurements on natural ferritin have shown that the 
cores-jdo iiave a small net magnetic moment (~100's of 
/i5)JIlE3"E^ presumably due to the preferential popula- 
tion of one magnetic sublattice during t.hi? formation of 
the particles. Recent theoretical workEi'tj predicts that 
an excess spin wlOO will have a small but appreciable ef- 
fect on the MQC frequency of an antiferromagnetic parti- 
cle. Motivated by this prediction and the opportunity to 
test long standing modclsO of the size dependence of the 
excess moment in small antiferromagnetic particles, we 



have measured the excess moments of several artificially 
synthesized ferritin samples. 

Occurring in a wide range of plants, animals, and bac- 
teria, ferritin consists of an organic hollow spherical shell 
with an inner core ~80 A in diameter. This protein ab- 
sorbs Fe ions through channels in its shell where they nu- 
cleatejnto an insulating antiferromagnetic crystal (T^r = 
240iiO) similar to ferrihydrite. Natural ferritin can con- 
tain at most 4500 Fe(III) ions, and typically contains an 
average of 2000. Using synthetic chemical techniques, it 
is possible to prepare samples in which each protein shell 
contains a fairly well-specified number of Fe(III) ions. 
The proteins used in these measixrenients have been ex- 
tensively characterized elsewhere.clEllla The samples have 
nominal loadings of n=100, 250, 500, 1000, 2000, and 
3000 Fe(III) ions (ionic moment fJ.Fe{iii) = 5.92/1^) per 
protein. Transmission electron microscopy (TEM) mea- 
surements of the mean and variance of the core diameters 
have been made for all but the smallest two samples, and 
published elsewhcrclj For magnetization measurements, 
a dilute (~0.5mg/ml) solution of a sample is dried on a 
polypropelene film and mounted on a twisted Cu wire 
in a commercial SQUID magnetometer. Measurements 
are made at temperatures ranging from T=4-300K and 
applied fields H=0-5 T. 

Typical M(H) curves for the n=2000 and n=3000 sam- 
ple are shown in Fig.l and reveal the presence of two com- 
ponents: one which saturates at large JI, and a second 
that is approximately linear in H. NeeO modeled small 
antiferromagnetic particles with net moments as an or- 
der parameter (essentially the Neel vector) possessing a 
magnetic moment as well as parallel and perpendicular 
susceptibilities X|| and x±- The energy of such a particle 
can be written as 
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where H is the applied field and ip is the angle between 
the Neel vector and H. The thermodynamic magneti- 
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zation M{H) = ksT-^ where Z is the partition 

function, can be calculated explicitly, but for H < ^i/x± 
(which holds for all npasurements here), we use the ap- 
proximate expression :E3l3 

M(H) = x\\H + (m + 2(x± - X\\)kBT/fi)L[^], (2) 

where L[x\ is the Langevin function. Because Eq.(||) 
assumes thermodynamic equilibrium, all of our mea- 
surements are made above the blocking temperature Tg 
(which is size dependentS) of each sample. This also jus- 
tifies our ignoring anisotropy energy terms in Eq. (|l|) . 

Figs. 1(a) and 1(c) show a fit to Eq.(||) for the data 
from the ?t.=2000 and n=3000 samples a few Kelvin above 
Tb- The fit fails in the region where M{H) has the most 
curvature, as there is a non-negligible spread in the par- 
ticle sizes in each sample. This results in a departure 
from Eq.(^ which is difficult to model without know- 
ing the precise size distribution in the sample. However, 
we can exploit the fact that the susceptibility near zero 
field is simply the average susceptibility of the sample. 
We first subtract the component of M(H) which is linear 
at high fields where the Langevin function has saturated 
to remove both the contribution from the antiferromag- 
netic bulk (the first term in Eq.(^) and any diamagnetic 
background. The data is then normalized so that the 
saturation moment is unity and fi is extracted from the 
slope of the low field data, n/ {3kBT). This approach has 
the advantage of being insensitive to the total number of 
particles in the sample as well as to the details of its size 
distribution since it returns the average fj, of the sample. 
Figs. 1(b) and 1(d) show fits to the high- and low- field 
data from Figs. 1(a) and 1(c) using Eq.(||). The fits rep- 
resent the theoretical M(H) of a particle with fj, equal to 
the average of the sample. 

This procedure can be used at any temperature above 
Tb (i.e., whenever the magnetic moments are in thermal 
equilibrium on the time scale of the measurement). We 
repeated measurements of the type shown in Fig.l over 
a temperature range from roughly Tb to 4Tb. In the 
samples with larger cores (Fig. 2(a)), we find a weak tem- 
perature dependence of the extracted /z, increasing by 
^--^20% over a factor of 4 in temperature. For the small- 
est cores (Fig. 2(b)), the temperature dependence of /i is 
quite strong, increasing by roughly a factor of two over 
the same range. A similar trend has b een , observed in 
natural ferritin as well as NiO particles .EHlEJ The reason 
for this temperature dependence is not clear. It should 
be noted that the model behind Eqs.(|l|) and (||) does 
not take into account any of the microscopic phenomena 
which might alter the properties of a small antiferromag- 
netic particle except in as much as these effects, can be 
modeled by /i, X|| X±- Weaker exchange Jl3 strong 
radial anisotropics, til and frustrationcj can exist at the 
surfaces of such particles, and may be responsible for the 
observed temperature dependence of /j,. Multiple sub- 
lattices can also exist in very small particles of a mate- 
rial which is antiferromagnetic in the bulk.cl At the low 



concentrations used here the typical interparticle dipolar 
fields should be well below 1 G, too small to account for 
these effects. For comparison with MQC, which is only 
observed below 200 mK, we extrapolate the linear tem- 
perature dependence shown in Fig. 2 to T=0 in order to 
extract the relevant excess moment. 

The result of this extrapolation is shown in Fig. 3, 
where ^ is plotted vs. n. The vertical error bars repre- 
sent the combined effects of the reproducibility between 
identically prepared samples and the uncertainties in the 
linear extrapolations of Fig. 2, '^10%. Because the ex- 
tracted /i corresponds to the excess moment averaged 
over the sample, the horizontal error bars do not rep- 
resent the variance in particle size, but rather the un- 
certainty in the mean particle size, estimated from the 
discrepancies between the nominal loading and the par- 
ticle size measured by TEM ('^20%). This is probably 
an underestimate in the case of the smallest two samples. 
The data are fit by a 0.56 ± 0.05 power law. If one plots 
fi in units of HFe(iii), the coefhcient of the power law 
fit is 1.15, quite close to unity. Neel has suggested three 
models of imperfect sublattice compensation in small an- 
tiferromagnetic particles.ii3 In the first, an antiferromag- 
netic particle has a surface consisting of sites belonging to 
one sublattice only, giving ^ — en's ^p^^^ jjjy The propor- 
tionality constant c is roughly 4 for the platonic solids. 
In the second model, the surface sites are distributed 
randomly between the two sublattices; then one has a 
random walk over the surface, and so = c'^n's ^p^(^jjjy 
For a particle surface of fractal dimension (as predicted 
by models in which the ferritin core is formed by diffu- 
sion limited aggregationtj) , the random walk over the 
surface can give any power law from 1/3 to 1/2. Lastly, 
if all the ions (as opposed to merely the surface ones) 
are randomly distributed between the sublattices, then 
fi = n2 This can occur, for example, if there is 

a non-stochiometric replacement of some magnetic ions 
with non magnetic ions. This prediction (which has no 
free parameters) is plotted in Fig. 3 as a dashed line. 
The agreement between this last prediction and the data 
might be strong evidence for the random population of 
the sublattices throughout the volume of the particles. 
We note, however, that for particles with 100 to 3000 
magnetic ions the discreteness of the lattice, combined 
with any surface roughness, means that a disproportion- 
ately large number of sites will be located on the surface. 
Thus, it is not possible to determine whether the sub- 
lattice non compensation is a volume effect or a surface 
effect. That the power law is clearly much less than 1 
is strong evidence that the excess moment does not re- 
sult from canting of the sublattices. Measurements made 
on fully loaded natural ferritin (n = 4500) and partially 
loaded natural ferritin|4fir^ 2000) are in good agreement 
with a 1/2 power law.liilT^ 

We can now compare the measured values of the ex- 
cess spin and the MQC resonance frequency vmqc a-s 
a functions of n. An antiferromagnet strongly coupled 
to an uncompensated moment can be described by the 



2 



effective actionl!'^'^^ 



->E 



V 



H 



? sin^ I 



Ky sin^ 6 sin^ . 



+ cos^ \ +ihS dr 0(1 - cos 0), (3) 



where 9 and are the spherical coordinates of the Neel 
vector, V is the volume of the grain, > Ky > 
its magnetic anisotropics, 7 = and S is the 

magnitude of the excess spin. Let us define S'afm = 
hV ^y2KyX±/ fJ-B, the instanton action one would obtain 
for an antiferromagnet without an uncompensated mo- 
ment. 

We will use instanton techniqufts-|to calculate the tun- 
nel splitting. It has been showrJla that in the regime 
Ky <^ Kz instanton solutions have an approximate fre- 
quency 
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The tunnel splitting is then given by 
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and the crossover temperature to the quantum regime by 



ksT* = KyVh/Spciri- 



(9) 



We suppose Xfr-^ independent of loading and use the 
simple estimateOa x_l = IJ-%N/kBTNV. The volume 
of the magnetic core of a fully loaded grain is roughly 
that of a sphere with a diameter of 7.5 nm, and contains 
4500 ions. Tjy = 240 K in ferritin, which gives x± ~ 
5 • 10^^ emu/Gem^. 

The crossover termperature to the quantum regime has 
only been measuredQ in fully loaded grains, where T* = 
200 mK. In the absence of experimental data we will 
assume T* is also independent of loading. 

Armed with these values for x± and T*, we are left 
with two unknowns at each loading, the anisotropics Ky 



and Kz- We can deduce their values from Eqs. (||) and 
, using the experimental values for /z and vmqc listed 
in Table |. We see that the hypothesis that Ky <^ Kz is 
at least self-consistent, and that the anisotropics are not 
strongly size dependent. The only exception is the value 
obtained for Kz in fully loaded grains. We note that the 
measurement of /i in fully loaded natural ferritinE3 lies 
somewhat below the power law shown in Fig. 3, possi- 
bly as the result of slightly different sample preparation, 
interparticle effects (the samples were not diluted), or 
differences in the surface of a grain completely filling the 
spherical protein shell. The value of Ky increases some- 
what for smaller particles, consistent with, the trend in 
blocking temperature as a function of n,t3 though it is 
somewhat smaller than typical for antiferromagnetic par- 
ticles. We are not aware of any other measurements of 
the transverse anisotropy Kz in such systems. 

In conclusion, we have measured the excess moment of 
the antiferromagnetic protein ferritin as a function of the 
number of magnetic ions per protein. Using diluted sam- 
ples to ensure the absence of interparticle interactions, 
we find an approximately square root dependence of the 
excess moment upon particle size, in agreement with a 
simple model (which has no free parameters) of imperfect 
sublattice compensation. We use this result to compare 
recent theoretical work on the effect of an excess moment 
to earlier MQC and blocking temperature measurements 
in the same samples. 
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FIG. 1. Magnetization vs. applied field above the blocking 
temperature for samples with particle size n = 2000 Fe(lll) 
ions (a,b) and n = 3000 (c,d). The solid line is a fit to the form 
of Eq.(^, where in (b) and (d), only the low- and high-field 
data are fit, as described in the text. The insets in (b) and 
(d) are magnifications of the low-field data and fit. 
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FIG. 2. The excess moment fi (in units of fiB and the ionic 
moment yiFe(iii)) oi ferritin cores each with n Fe(III) ions as 
a function of temperature from Tb to roughly ATb for each 
sample. The dashed lines are linear fits to, in (a) the larger 
cores, and in (b) the smaller cores. The temperature scale is 
the same in both a) and b). 
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FIG. 3. The T=0 excess moment in units of /is as a func- 
tion of particle size n. The solid line is a power law fit giving 
an exponent of 0.56. The dashed line is the no free parameter 
prediction of the \i — iJ,Fe{iii) model described in the text. 
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TABLE I. Excess moment /i, MQC resonance frequency 
vmqc, and anisotropy energies Ky and for ferritin cores 



containing n Fe(III) ions. 
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